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Abstract: We examine whether the ATLAS detector has sensitivity to extra- 
dimensional scalars (as opposed to components of higher-dimensional tensors which 
look like 4D scalars), in scenarios having the extra-dimensional Planck scale in the 
TeV range and n > 2 nonwarped extra dimensions. Such scalars appear as partners 
of the graviton in virtually all higher- dimensional supersymmetric theories. Using 
the scalar's lowest- dimensional effective couplings to quarks and gluons, we compute 
the rate for the production of a hard jet together with missing energy. We find a 
nontrivial range of graviscalar couplings to which ATLAS could be sensitive, with 
experiments being more sensitive to couplings to gluons than to quarks. Graviscalar 
emission increases the missing-energy signal by adding to graviton production, and 
so complicates the inference of the extra-dimensional Planck scale from an observed 
rate. Because graviscalar differential cross sections resemble those for gravitons, it 
is unlikely that these can be experimentally distinguished from one another should 
a missing energy signal be observed. 
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1. Introduction 

The realization that gravitational physics could be as low as the TeV scale [[I], fj, 
has sparked considerable interest in the possibility of detecting extra-dimensional 
gravitons in upcoming accelerator experiments @, |j| There are now several types 
of models within this class, which largely differ in the masses which they predict 
for extra-dimensional states. In the extreme ADD scenario these Kaluza-Klein 
masses can be as small as rriKK ~ 10~ 3 eV, while in the warped RS picture 0, 
they are as large as the TeV scale. 

Although the graviton has grabbed most of the phenomenological attention, 
models with low string scales often also predict many other kinds of extra-dimensional 
physics which might also be detectable. For instance, many of the best-motivated 
models are inspired by string theory, and so often predict the gravitational physics of 
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the bulk to be supersymmetric and thus contain entire extra-dimensional gravity su- 
permultiplets. Indeed the supersymmetry of the gravitational sectors of these models 
may be one of their great strengths, perhaps helping to control the contributions of 
quantum corrections to the effective 4D cosmological constant 0. 

Since the relevant supersymmetry is extra-dimensional, from the four- dimensional 
perspective the graviton supermultiplet forms a representation of extended super- 
gravity (i.e. supergravity with more than one supersymmetry generator). As such, 
the graviton multiplet in these theories can contain several spin-3/2 gravitini, spin-1 
gauge bosons, spin-1/2 fermions and scalars, in addition to the graviton. For example 
a typical gravity supermultiplet in 6 dimensions already contains 4D particles having 
all spins from through 2. Although much less is known about the phenomenological 
consequences of these other modes — see, however, || - generically they may be 
expected to have masses and couplings which are similar to the graviton's since they 
are related to it by extra-dimensional supersymmetry. 

In this paper we focus on the properties of the graviscalar, which we take to 
mean an honest-to-God extra-dimensional scalar which is related to the graviton 
by supersymmetry. (We do not, for example, mean a 4D scalar which arises as an 
extra- dimensional component of the metric tensor itself — although many of our 
results will also apply to such a particle [fL0| .) Our goal is to identify the relevant 



couplings of such a scalar, and to use these to explore its experimental implications. 
In particular, we identify how the production of such a scalar can compete with the 
predictions for graviton production. 

For phenomenological purposes we concentrate on the implications of real gravis- 
calar production, rather than calculating the consequences of its virtual exchange. 
We do so for the same reasons as also apply to virtual exchanges of gravitons || 
and graviphotons ||: their virtual exchange cannot be distinguished from the effects 
of local interactions produced by other kinds of high-energy physics (such as the 
exchange of massive string modes) ||. 

We organize the presentation of our results as follows. In the next section we de- 
scribe the most general possible low-energy couplings of an extra- dimensional scalar 
with quarks and gluons, and summarize the domain of validity of an effective-field- 
theory analysis. Section 3 then uses these couplings to compute the relevant cross 
sections at the parton level, and for jet plus missing energy production in proton- 
proton collisions. Section 4 summarizes the results of numerical simulations based on 
the cross section of section 3. It compares the predicted production rate both with 
the expected Standard Model backgrounds, and with previously-calculated graviton 
emission rates. Our conclusions are then briefly summarized in section 5. 

2. Low-Energy Graviscalar Couplings 

The class of models of present interest are those for which all Standard-Model par- 
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tides reside on a 4- dimensional surface (or 3-brane) sitting in a higher- dimensional 
space. We would like to know how such particles can couple to higher- dimensional 
scalars which are related to the graviton by supersymmetry. In particular, since we 
look for direct experimental signatures in accelerators, we concentrate on trilinear 
interactions involving a single higher-dimensional scalar and two Standard-Model 
particles. 

At low energies we may parameterize the most general such couplings by writing 
down the most general lowest-dimension interactions which are consistent with the 
assumed low-energy particle content and symmetries of the theory. In this paper, 
we are concerned with the production of scalars in parton interactions. The most 
general couplings to quarks and gluons are given by the following Lagrangian: 

£eff = d M (f>(x, y) d M (p(x, y) 

-5 n (y) Y,%( x )(9 + Wsls)ii%{x)<t>{x) (2.1) 
- Q 

+ c g G^{x)Gl v {x)4>{x) + b g e^Gl l/ (x)G a Xp (x)c/ ) (x) 
+ c 1 F^{x)F„ v {x)(t>{x) + b 1 e» l ' Xp F l , 1/ (x)F Xp (x)<f)(x) 

with arbitrary dimensionful coupling parameters g^, (g 5 )ij, c g , b g , c 7 and b 1 . 
The indices i, j = 1, 2, 3 here label the Standard Model's three generations. 

In what follows it is convenient to define dimensionless couplings by scaling out 
the appropriate power of the reduced Planck mass in D dimensions, 1 according to 
9ij = ~9ijM r p 2 , c g = CgM 1 ^ 1 ^ 2 and so on. It should be kept in mind when doing 
so, however, that the fermion interactions are not SUl(2) x Uy{l) invariant and 
so the size to be expected for the couplings g^ and (g 5 )ij depends strongly on the 
way in which the electroweak gauge group is broken in the underlying theory which 
produces them. In particular, if the new physics is not involved in electroweak 
symmetry breaking then there is a natural suppression of the fermionic dimensionless 
couplings, gij ~ v/Mp, where v = 246 GeV. (If the relevant dimensionless couplings 
of the underlying model are similar in size to Standard Model yukawa couplings, 
then this suppression can be even smaller, although they need not be this small in 
all models.) In principle, if the extra-dimensional physics were itself the sector which 
broke the electroweak gauge group, even the suppression by powers of v/ Mo might 
not be present, although in this case the new-physics scale cannot be very large 
compared to v 

In these expressions, the coordinates x M describe the 4 dimensions parallel to the 
Standard-Model brane and y m similarly describes the n transverse dimensions. The 



1 More precisely, M D is the reduced Planck mass in D = 4 + n dimensions, defined in terms of 

2 — D 

the D-dimensional Newton's constant by 87r Go = M D 
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brane itself is assumed to be located at the position y m = in these extra dimensions. 
F^y denotes the usual electromagnetic field strength, and G°L, is the non-abelian gluon 
field strength, G° v = d ll Ql-d v G*+g z f abc G b ll Gl, with g 3 denoting the QCD coupling 
constant. Finally, \I/q generically denotes any of the Standard- Model fermion mass 
eigenstates, whose electric charge is denoted by Q. 

Several comments concerning these effective interactions bear emphasis: 

• As written, the fermion interaction of eq. ( [2.1| ) appears not to be invariant under 
electroweak gauge transformations. This is because we have already replaced an 
explicit factor of the Standard Model Higgs field with its vacuum expectation value, 
v = 246 GeV, and rotated to a fermion mass eigen-basis. Although the resulting 
couplings, gij and (gs)ij, can in principle be off-diagonal and so involve flavour- 
changing neutral currents, we will assume these not to arise since they would be 
strongly constrained if present. 

• Keeping in mind that in real models the scalar in question typically lies in a su- 
permultiplet with the graviton and so couples with similar strength, the interactions 
written above can be regarded as of leading order in inverse powers of M^, and 
are expected to break down once physical energies approach this scale. One excep- 
tion would be the possibility of a direct mixing term between the extra- dimensional 
scalar and the Standard Model Higgs. Since the following analysis concentrates on 
the phenomenology of jets plus the graviscalar, we ignore this kind of 0-Higgs mixing 
in what follows. For similar reasons we also ignore graviscalar couplings to the Z 
and W bosons. 

• Derivative couplings to the fermions are not included here since such terms are not 
independent, as they can be re- written in the form given above by performing a field 
redefinition and are expected to be small |12|| . 

• The factor 5 n (y) in eq. ( |2.2| ) expresses explicitly the broken translation invariance of 
the bulk due to the presence of the brane. Consequently momentum transverse to the 
brane is not conserved, allowing bulk particles to be emitted into the extra dimensions 
even if all of the initial particles of the interaction were themselves confined to the 
brane. Physically, the unbalanced transverse momentum is absorbed by the recoil of 
the brane itself, with no energy cost because of the brane's enormous mass. 

• Our effective Lagrangian, as written in eq. (|2~l| ), applies equally well to both 
warped |§ and unwarped || models. The spectrum of effective 4D masses emerges 
from the expansion of the bulk scalar field in terms of a basis of Kaluza-Klein 
(KK) modes: (fi(x,y) = J2k ( Pk( x ) u k(y)- Below, we assume large extra dimensions, 
allowing us to integrate over the phase space of all bulk particles using a flat geom- 
etry (for a more explicit computation of the limitations of this approximation, see 
ref. 0.) 
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3. Physical Predictions 



The effective Lagrangian of eq. Q2.1|) can be used to compute the cross-section at 



tree-level for production of graviscalars in p-p collisions at the Large Hadron Collider 
(LHC). We shall restrict our attention to hadronic production mechanisms since 
these are expected to dominate at LHC energies, and so might be expected to give 
the clearest signal. For these purposes our interest at the parton level is therefore in 
the reactions qq — > <pX, gg —>■ <pX and qg —>■ <f)X, where X represents a Standard- 
Model final state such as a single well-defined jet. The physical signal which such a 
reaction would produce is a well defined jet plus missing energy as the graviscalar 
escapes into the extra dimensions. 

3.1 Parton-Level Cross-Sections 

Using these Feynman rules (given in Appendix), we next compute the parton-level 
cross section for producing a hard quark or gluon plus a graviscalar in the final state. 
There are three processes which are relevant: qq — > g(p, qg — > q<j) and gg —>■ g<p. 

It is convenient to divide the higher- dimensional graviscalar momentum vector, 
^4+ n with M — 0, ...,3 + n, into its continuous 4-dimensional components, l\ with 
/i = 0, 3, and its quantized n-dimensional components, L m with m = 4, 3 + n. 
The total squared- momentum for a massless bulk graviscalar is then £ 2 +n = i\ + L 2 = 
0, which we write as £ 2 = —M 2 , where M is the effective 4-dimensional mass due 
to the particle's motion in the extra dimensions. As mentioned earlier, we perform 
our calculations in the approximation where the length scales associated with the 
extra dimensions are much larger than the wavelengths of the partons involved, since 
this is a good approximation for the applications of interest and allows a tractable 
treatment of the graviscalar phase space in the extra dimensions. (However, as a 
consequence, our further results do not generalize to warped models.) With the 
approximation the quantization of L m is not important, and sums over this variable 
may be approximated as integrals. The corresponding integration measure is then: 



After integration over the angular degrees of freedom, Q n , the final phase-space 
measure used for the graviscalar bulk momentum becomes: 



Evaluating the Feynman graphs of Fig. (^) (see Appendix) to obtain the parton- 
level differential cross sections gives the following results. The quark-annihilation 
channel is obtained by evaluating graphs (a) through (c), which give: 






(3.2) 
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da(qq 



a s (2n) n ' 2 {M 



2\(n-2)/2 



dtdudM 2 



18r(f)MSs 2 



(g 2 + g!) f2M 2 s + (u + t) 



(2 7r )2n/(2+n) 



Ml 



4(c 2 + 6 2 ) ft 2 + u 



ut 

\2\ 1 



(3.3) 



5(s + i + u + M 2 ) . 



Notice that helicity conservation precludes graph (c) interfering with graphs (a) and 
(b) in the limit where parton masses are neglected (as we assume). 

The quark-gluon scattering contribution is similarly obtained by evaluating the 
graphs (d), (e) and (f), to give 



da{qg -> q(f>) 
dtdudM 2 



a s (2vr)"/ 2 (M 2 )(™- 2 )/ 2 



48T(f)M™ s 2 



(9 2 + 9l) fu 2 + M 



( 27r )2n/(2+n) 



St 



(3.4) 



4(c 2 



M 2 D 



g + s 2 
u 



5(s + i + u + M 2 ). 



As before, the neglect of parton masses ensures the non-interference of graph (f) with 
graphs (d) and (e). 

Evaluating the gluon fusion graphs, (g) through (j), gives 



da(gg -> g<j>) 3 a s (2n) n / 2 (M 2 )( n ~ 2 V 2 ( (c 2 + b 



dtdudM 2 



16T(f)s 3 tM 



Ml +2 



+12 si uM 2 + -{i 2 u 2 + i 2 s 2 + u 2 s 2 ) 
2 



u + i) 4 + ( u + s) A + (i+s) 4 

5(s + i + u + M 2 ). (3.5) 



Notice that both the expressions for the gluon-fusion and quark-annihilation pro- 
cesses are invariant under the exchange t <-> u, as is expected on general grounds 
from the charge-conjugation invariance. The mass scale Mu used here is related to 
the scale M D defined earlier by M™ +2 = (2vr) n M^ +2 . 

Since these expressions depend only on the dimensionless coupling combinations 
g 2 + g\ and c 2 + b 2 , in what follows we set g 5 = b = and choose c > and g > 
without loss of generality. 



3.2 Proton-Proton Cross Sections 

The cross section for proton-proton collisions is obtained from the parton-level results 
just calculated in the usual way, by convoluting with the parton distribution func- 
tions, fi(x,Q 2 ). For our later analysis we compute the cross section for the process 
pp — > + jet, which has the form: 



d Xl dx 2 dtdM 2 fi( Xl , Q 2 ) fj(x 2 , Q 2 ) da ^ 3 ~* 



dtdM 2 



(3.6) 



s=xi X2 E c 
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where the sum on i,j runs over the types of partons available in the proton and X 
corresponds to an energetic quark or gluon. We generate events randomly to which 
we assign weights by performing the integration. These events are then accepted 



or rejected according to their weight by the generator PYTHIA [14], which also 
performs hadronization. 

The phase-space integrals are performed subject to the following constraints. 

• We require the final jet's transverse momentum to satisfy P% > P 2 ut , where P cut 
(= E™jl t ), is a minimum value which we specify below. Using the kinematical 
relation P% —tu/s, where u = M 2 — i— s, this leads to the following integration 
limits for t: t min < i < t max , where 

trmn = \ [(M 2 - s) - yj (M 2 - - AP^S ] (3.7) 

and t max = l - [(M 2 -s) + \J '{M 2 - s) 2 - 4P 2 ut s ] . (3.8) 

• Energy-momentum conservation implies the following upper limit for M 2 

< M 2 < M 2 max = s - 2V§P cut . (3.9) 

• Finally, the integration limits on the parton energy fractions are: 

s ■ AP 2 

°min rut ^ ^ 1 in i r\\ 

x min = = — 2 ^<x 1 ,x 2 <l, 3.10 

s s 

where as usual s denotes the invariant initial energy of the proton-proton col- 
lision. 

We have checked our numerical integration by recomputing the differential cross- 
section for graviton production, using the parton-level cross sections computed by 
Giudice et al 0, and comparing with their results. 

3.3 Validity of the Low-Energy Approximation 

As discussed above, our expressions for the parton-level reaction cross sections are 
only valid for parton energies well below the cut-off scale, Mo- But since we can only 
specify the initial proton energies (and phase-space cut-offs like P cu t) , we cannot know 
for sure whether the numerical integration includes parton reactions which carry 
energies which are too large. If so, our calculation becomes unreliable in that part of 
phase space for which the probability of very-energetic parton processes cannot be 
neglected. In this section we define a proton-level criterion for estimating the extent 
to which high-energy parton processes pollute our calculations in various parts of 
phase space. Our goal in so doing is to be able to choose upper limits for quantities 
like P cu t which minimize this pollution. 

In order to do so, we compare in Fig. ([I]) the total proton-proton cross-section, 



a, calculated in the following two ways HL ITB| : 
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1. We calculate using the above parton cross sections for all possible parton en- 
ergies. 

2. We calculate using the above parton cross sections only if \/| < M D , and set 
the parton-level cross-section to if vl > Mb- 

The bottom panel of Fig. (|I|) shows this comparison when only the quark-graviscalar 
couplings, g,gs, are nonzero. The top panel shows the result assuming only the 
gluon-graviscalar couplings, c, b, do not vanish. The various curves plot a against 
E™ l jl t = P cut for different choices for s max (with the effective dimensionful couplings 
~g and c held fixed at some arbitrary reference value) When the two curves start 
to deviate, high-energy parton contribution are significant and we do not trust our 
calculation. 

We use these curves to define the maximum value of P cu t which we may trust, 
given a value for M D (and so also for the effective graviscalar couplings). Quantita- 
tively, we fix P cut by demanding that the curves not differ by more than 10%. As is 
clear from the figure, the value of P cu t which is obtained in this way is smaller for 
the gluon-graviscalar couplings (g = g^ = 0) than for the quark-graviscalar couplings 
(c = b = 0), and we use the lower of the two in the following calculations. 

Notice that these two plots also indicate that for numerically equal couplings, 
the gluon-graviscalar couplings dominate the cross-section at high energy, while the 
low-energy regime gets bigger contributions from the quark-graviscalar interactions. 
This is as expected from the effective Lagrangian, since the gluon terms involve an 
addition derivative (and so an additional power of E/Mb in cross sections) relative 
to the quark terms. 

For fixed P cu t, we may ask how large the effective graviscalar couplings may be 
without introducing more than a 10% error into our calculations. This is illustrated 
in Figs. (0), which show the smallest value of Mb which is permitted given a choice 
for E™j1t = Pcut- This lower bound on Mb amounts to choosing an upper bound for 
the effective dimensional couplings, g and c. 

In what follows we calculate the proton-proton cross sections assuming a value 
for Pcut, and so the condition Mb > M D VLn determines the domain of validity of 
our calculations. As discussed above, we determine M% m using gluon-graviscalar 
couplings, since this is the stronger requirement. From this we find the minimum 
values for ~g and c for which a graviscalar signal would be observable above the 
statistical Standard Model background at the 5 o level. 

4. Simulations 

As mentioned above, the processes discussed in Sect. ?? were implemented using 
PYTHIA as external processes. Parton flavors were properly assigned in each event 
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Figure 1: Total jet + nothing cross-section for graviscalar production at LHC to evaluate 
the validity limit of the phase space volume when: top) only the graviscalar-gluons cou- 
pling is present; bottom) only the graviscalar-fermion coupling is present. The curves are 
normalized to the cross section for graviton production at a value of E™™ t of 500 GeV. 



according to the CTEQ 5L parton distribution functions evaluated at the renormal- 
ization scale Q 2 = | M 2 + p\ and the color flow between those parton was applied. 
ATLAS detector effects were incorporated using the fast Monte Carlo program ATL- 



FAST . 

For the purpose of comparing with Standard Model backgrounds, we must choose 
reference values for the couplings g and c as well as for the number of extra dimensions 
n. For these purposes a useful choice is a set of couplings for which the low-energy 
approximation applies and for which the proton-proton cross section is the same size 
as the cross-section which was determined as being what is required for a 5 a discovery 
for graviton production [[ll|. The reference dimensionless effective couplings, when 
Pcut = 500 GeV, Md = 5 TeV and for n = 2 extra dimensions, are then g ~ 0.70 
or c ~ 0.41. These couplings are reasonable from the point of view of the effective 
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Figure 2: Evaluation of the minimum (4+n)-dimensional Planck scale for which the 
effective model is valid, at the 90% level, given a cut E^ l ^ eV as explained in the text 



theory because the g and c we obtain in this way are smaller than unity. Notice 
that the existence of such couplings shows that there exist scenarios for which the 
graviscalar production would be as important as graviton production. 

We remark that the following analysis is independent of this choice of reference 
values. 

4.1 Standard Model Backgrounds 

If graviscalars are produced in association with a jet in proton collisions the event 
can be found by searching for a jet plus missing energy: pp — > jet+^ T . The Standard 
Model background to this process arises from events having neutrinos in the final 
state. The principal backgrounds of this type and their cross-section in the phase 
space regions E^ n et > 500 GeV and E™™ t > 1000 GeV are given, following ref. 



in Table ([[]). For comparison, the graviscalar production cross-section using our 
reference couplings — 2 extra dimensions, Md = 5 TeV, E™ l jl t = 500 GeV and 
(c, g) = (0.41, 0.70) — is a = 156 fb. 
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Processes 


cross-section (fb) 


500 GeV 


1000 GeV 


pp — >jet+Z(— > vv) 


278 


6.21 


pp ^jet+W (— > ev e ) 


364 


8.57 


pp ->jet+VK(— > \iv^) 


364 


8.51 


pp — >jet+W(-^ rv T ) 


363 


8.50 



Table 1: S.M. background to the graviscalar production at ATLAS and their cross-section 
for different phase space volume. 

4.2 Analysis 

We now determine a set of experimental cuts with which we can find the minimum 
values of g and c for which a 5a discovery is possible. Jets are reconstructed using 
the cone algorithm with a cone radius AR = 0.4. In the detector ATLAS, leptons are 
detected if they are emitted in the range of pseudorapidity —2.5 < r\ < 2.5. Leptons 
are defined as isolated if the energy deposited by other particles in a cone of radius 
AR = v/(A0) 2 + (At/) 2 is less that 10 GeV. We impose the following two cuts. 
As a first cut we require: 

Cut 1: No isolated lepton (electron or muon) with p T > 6 GeV is allowed 
in the event. 

This eliminates most of the W — > ev e and W — > fiv^ events, leaving only those 
for which the leptons are not properly reconstructed in the detector. It does not 
eliminate events like W — > tv t in which the r decays hadronically, however in this 
case we expect to also have an energetic low-multiplicity jet which is opposite, in the 
azimuthal plane, to the principal jet. We see in Fig. (|3p that a cut on the difference 
in azimuthal angle between the two most energetic jets can eliminate a significant 
fraction of this r background. 

We therefore impose the second cut: 

Cut 2: We keep only events for which \ip^ — (fj 2 \ < 2.285 radians. 

This cut is chosen to maximize the significance of the remaining signal. 

Fig. (|) shows the relative contribution of each process to the total background. 
After the above cuts, the most important background is pp — > ]etZ — > jet vv. Ta- 
ble (HI) makes this more explicit, by breaking down the background and comparing 
it and the signal before and after the cuts are applied, assuming an integrated lumi- 
nosity of 100 fb" 1 . For this integrated luminosity — which corresponds to one year's 
running at the nominal LHC luminosity — we therefore expect a total of 36700 
background events to remain after cuts. 
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Figure 3: Distribution of the difference in the azimuthal angle between the two most 
energetic jets of an events for: (top) each of the background; (bottom) the signal and the 
total background on top of it. 



Is this good enough? We estimate the number of signal events required for a 5<r 
discovery using the following significance criterion: 

S >5, (4.1) 



VS + B 



where S and B are respectively the number of signal and background events. For 
this many background events we therefore have a 5cr discovery if more than 970 
graviscalar events are detected, i.e. if the total cross-section for the process is: 
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Processes 


E™ in > 500GeV 


# events after cut 1 


# events after cut 2 


jet+Z(^ vv) 


27760 


27100 


24940 


jet+W(-> ev e ) 


36420 


5224 


1430 


jet+W(-> nv„) 


36370 


957 


866 


jet+W(-» rv T ) 


36330 


24600 


9459 


jet+Graviscalar 


30960 


30090 


27720 



Table 2: Number of signal (Md = 5 TeV, n = 2, g = 0.70 and c = 0.41) and background 
events that survive each cut for an integrated luminosity of 100 fb . The cuts are defined 
in the text. 




600 800 1000 1200 1400 1600 1800 



Figure 4: Contributions of different processes to the total background after application 
of the cuts, for an integrated luminosity of 100 fb -1 . 

cr(pp — > jet + 0) > 10.9 fb. For n = 2 this corresponds to the effective couplings: 

g> 0.18 TeV" 1 if c = 0, (4.2) 

or 

c > 3.2 x 10" 1 TeV" 2 if g = . (4.3) 

Fig. (P) plots the cross section as a function of missing energy, for two choices 
of couplings. The choice in the second panel corresponds to the 5a discovery limit, 
and shows that the discovery would be due to an excess of events in the distribution 
of missing transverse energy at high energies. 
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600 800 '000 1200 1400 1600 1BO0 




600 800 '000 1200 1 400 1600 1800 



E T mlss 

Figure 5: Distribution of missing Et for background and signal when: (top) the cross- 
section for graviscalar production is the same as the graviton one; (bottom) this cross- 
section is at its discovery limit. 



Notice that the cross section required for discovery is determined by the back- 
ground rate and so is the same for all choices for n, the number of extra-dimensions. 
In fact, as can be seen in eq. fl3.5|) , the t and u dependence of the cross-section does 
not depend on the number of extra dimensions. Only the effective coupling constants 
and the energy dependence depend on n, through the overall factor (M 2 )*™ -2 )/ 2 . Since 
the angular distribution depends only on t and u, Cut 2 has the same effect for all 
possible n (as does Cut 1). 
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4.3 Results 

We now turn to the central question: What range of effective couplings are likely to 
be detectable at the LHC? 

We have seen that any determination of the reach of LHC must be made relative 
to a choice for E™ l jl t = P cu t, since this plays a role in the reliability of the entire 
theoretical calculation. From Fig. (Q) we see that the choice EJp % J l et = 500 GeV implies 
that the cross-section is sensitive to high-energy parton processes at less than the 
10% level, provided M D > M% in , where M% in = 3.60, 4.30, 4.85 and 5.70 TeV for 
n = 2, 3, 4 and 6 extra dimensions respectively. Given the value for M™ 71 we then 
determine what values of couplings produce an observably large cross section (i.e. 
a > 10.9 fb). 

Suppose (#o6g,c^) are a pair of dimensionful couplings which each by itself 
produces a 5a signal (in n extra-dimensions) above the Standard Model background. 
The LHC then can detect couplings which lie in the intervals 

l>9>9^(M^ m ) n/2 if c = 
or 

i>c>4!(C) 1+n/2 if g = o. 

The upper limit in these inequalities expresses the theoretical criterion that the 
calculation only makes sense below the cut-off scale M™ n . 

Fig. (^|) shows the couplings which are accessible if both c and g are simulta- 
neously nonzero. The shaded regions indicate the range of couplings which are too 
small to have detectable effects for various choices for the dimension n. The poten- 
tially observationally-interesting couplings are those which lie outside the ellipses, 
but inside the box defined by c < 1 and g < 1. (Recall that since only c 2 and g 2 
enter into the cross sections, these plots should be interpreted as constraints on |c| 
and \g\.) 

We see that whether useful constraints are possible depends on the number of 
dimensions n, and fewer dimensions gives better reach. If n = 2 couplings outside the 
innermost region are potentially detectable. For n = 3 detectable couplings must lie 
outside the next-to-innermost region. Detection for the case n=4 is unlikely for the 
dimensionless constant g in the limit of small c, but is possible if c is also nonzero. For 
more then 6 extra dimensions, detection of any signal is unlikely since the minimum 
value of couplings needed for a discovery when n = 6 is (c, g) = (0.9,3.3). Any 
coupling in this region is far enough from the limits of validity of the calculation to 
have confidence in the result. 

An alternative way of expressing the potential ATLAS reach for a graviscalar 
signal is in terms of the value of the fundamental Planck scale to which the detector 
might be sensitive. It is bounded on the low side by the requirement that the effective 
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graviscalar-fermions coupling (g) 



Figure 6: Parameter region for graviscalar theory at second order in E/Md that allow 
testable and valid physical prediction at ATLAS, for different number of extra-dimension. 
The graviscalar-fermion and graviscalar-gluon dimensionless couplings are effectively com- 
binations (g 2 + gl) 1 / 2 and (c 2 + c|) 1//2 



ndim 


M min ( Te y^ 


M rnax (r^y) 


2 


3.2 


14.00 


3 


3.8 


6.45 


4 


4.4 


1.45 



Table 3: Sensitivity of the ATLAS detector to the fundamental Planck scale Mo through 
the discovery of a graviscalar signal, for c = and for an integrated luminosity of 100 fb _1 . 
For n > 4, observation of a signal is not possible. 



theory be a good approximation, and on the high side by the condition that the signal 
be detectable. Denning the effective upper limit by 

M% ax = (2vr)^ x min 

there are prospects for detection when M™ n < Mp ax . Table (|3|) summarizes these 
results, with M^ ax calculated using the worst case: c = 0. 

In the more optimistic limit g = 0, the maximum value of the fundamental 
Planck mass to which the ATLAS detector is sensitive increases from 6.45 to 9.50 
TeV when n — 3, and from 1.45 to 7.55 when n — 4. These results are summarized 
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-2/(n+2) 



(4.4) 
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in Table (Q). We again find that n = 6 is the limiting case since for g = 0, we have 
Mr = 5.8 TeV « MF™ n = 5.7 TeV. 



ndim 


M min 


M rnax ( Te y) 


2 


3.60 


14.10 


3 


4.30 


9.50 


4 


4.85 


7.55 


6 


5.70 


5.80 



Table 4: Sensitivity of the ATLAS detector to the fundamental Planck scale Md through 
the discovery of a graviscalar signal, for g = and for an integrated luminosity of 100 fb _1 . 
For n < 6, observation of a signal is possible. 



A comparison with the results obtained from graviton emission [|IBJ is also in- 
structive, although some care must be taken in so doing because the graviton results 
were obtained using a more restrictive phase-space cut {E™ m > 1 TeV), a different 
criterion for defining the validity region of the model and with a more conservative 
statistical estimator (S > V7B). Tables (§) and (|) compare the sensitivity of AT- 
LAS to Md as computed using graviscalar and graviton production, using these more 
conservative criteria. The two tables differ in their choice of either g = or c = 0. 

Table (|5]) also shows the existing non-accelerator limit on Mp, taken from ref. |L7 
(see also ||). Unlike the situation for gravitons (which couple universally) these 
astrophysical bounds are more model-dependent when applied to graviscalars. This 
model dependence arises because they directly bound the couplings of KK modes 
to electrons and photons, and so need not directly apply to the gluon and quark 
couplings of most interest for colliders. 



c=0 


Graviton 


Graviscalar 


limit from cosmology 






M max 


M min 


M max 


M"g in (A) 


M% in (B) 


n=2 


~4.0 TeV 


7.5 TeV 


4.35 TeV 


5.45 TeV 


O(90) TeV 


~ 10 TeV 


n=3 


~4.5 


5.9 


4.85 


3.65 


5.0 


0.8 


n=4 


-5.0 


5.3 


5.35 


3.20 


< 4 


<1 



Table 5: With c=0: comparison of the sensitivity of ATLAS to Md for graviscalar and 
graviton signals under the conditions E™ in > ITeV, -J= and with indirect constaints from 
cosmology. The integrated luminosity is 100 fb . For n > 3, observation of a graviscalar 
signal is not possible since M™ n > M^ ax . For the cosmology bounds, Scenario A means 
limits to neutron star heating by KK-decays, while scenario B corresponds to bounds from 
the cooling of SN1987A by KK-mode emission. 



We see again that the scenario where g — > gives the best case for detection, and 
this is competitive with the graviton result. We also see that although accelerator 



g=0 


Graviton 


Graviscalar 




M min 


M raax 


M min 


M max 


n=2 


~4.0 TeV 


7.5 TeV 


4.65 TeV 


10.20 TeV 


n=3 


~4.5 


5.9 


5.15 


7.75 


n=4 


-5.0 


5.3 


5.60 


6.50 



Table 6: With g=0: comparison of the sensitivity of ATLAS to Md for graviscalar and 
graviton signals under the consitions E™ m > ITeV, and with indirect constaints from 
cosmology. The integrated luminosity is 100 fb _1 . For n < 4, observation of a graviscalar 
signal is possible. 

experiments are most sensitive to lower n, for quark couplings these may be pre- 
empted by the non-accelerator bounds. 

The difference between the cases c = and g = show that ATLAS is likely 
to be only weakly sensitive to the graviscalar Yukawa couplings (especially keeping 
in mind these are naturally expected to be at most of order v/Mp, as explained in 
section 2), and a discovery is more likely to come from gluon-graviscalar couplings. 
However, once a signal is seen we are unlikely to be able to decide directly on the 
relative importance between g and c. Therefore, even if the discovery of a signifi- 
cant graviscalar signal at ATLAS should turn out to be possible, it is unlikely to 
completely fix its couplings. 

4.4 Graviton-Graviscalar Confusion 

Should a missing-energy signal be seen at the LHC, how does one tell if it is due to 
gravitons or graviscalars? We do not yet see a way to do so, despite the difference 
in their spin, for the following reasons. 

• The graviscalar production cross-section has an energy dependence which is 
similar to the graviton one, precluding the use of PTj e t,E r or an y °ther function 
of energy to discriminate the two. 

• Parton-level discriminants are not likely to be of practical use, because the 
center of mass energy of the hard scattering is not known in a pp collider such 
as the LHC. Furthermore, the final state we consider consists of a single jet and 
missing transverse energy, so it is not possible to reconstruct the longitudinal 
component of momentum of the system of interacting partons, nor their angu- 
lar distribution in the center of mass, nor their forward-backward asymmetry. 
Even if this were possible, we have checked that the discrimination between 
the shapes of the graviscalar and graviton differential cross sections is difficult 
even at the purely theoretical parton level. Only gluon fusion processes lead 
to a small difference. 
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5. Conclusions 



We have computed the rate for the production of extra-dimensional scalars (as op- 
posed to components of extra-dimensional tensors which look at low energies like 4D 
scalars) in pp collisions. Such particles arise in virtually all supersymmetric higher- 
dimensional theories, and our work complements previous studies of gravitons [|], 16 1 



and of extra-dimensional vectors |§. Because of the way we compute our phase space 
integrals, our study applies to large-extra-dimensional (ADD-type) models and not 
to warped (RS-type) models. 

We find that the cross sections for the reaction pp — > <fi + jet are similar in size 
and shape to those for graviton production, although the competing non-accelerator 
constraints on the couplings can differ because graviscalars need not couple univer- 
sally (unlike gravitons). We used simulation codes tailored to the ATLAS detector, 
and conclude that ATLAS can be sensitive to graviscalar couplings, provided there 
are less than 6 extra dimensions. The sensitivity improves with fewer dimensions, al- 
though so does the restrictiveness of non-accelerator bounds on the extra-dimensional 
Planck scale. Nontrivial windows of opportunity can be consistent with all bounds. 

Generically, pp collisions are more sensitive to graviscalar couplings to gluons 
than they are to couplings to quarks. Both couplings are turned on in our analysis, 
and we find that observable quark couplings often push the limits of validity of the 
effective- field-theory description. 
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7. Appendix 



To proceed the cross sections, we require the Feynman rules for the qq<f) } ggcj) and 
ggg<j) vertices which follow from the effective lagrangian of eq. ( j2.1[ ), which are: 



Gs 




Gs 
l 1 




= Ai[c(p.q)g^ - cp^q u - bt l , uaP p a q f3 ]5 ah 




= ^f abc [cg^{p P - q P ) + cg w {k v - p v ) + cg^q^ - fc M ) - be a ^ up (p a + q a + k a )} 
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Figure 7: The parton-level Feynman graphs which contribute to graviscalar production 
with an associated jet in proton-proton scattering. 



